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ABSTRACT 

We have carried out mid-infrared spectroscopy of seven Galactic proto- 
planetary nebulae (PPNs) using the Spitzer Space Telescope. They were observed 
from 10—36 /xm at relatively high spectral resolution, RRi600. The sample was 
chosen because they all gave some evidence in the visible of a carbon-rich chem- 
istry. All seven of the sources show the broad, unidentified 21 /im emission 
feature; three of them are new detections (IRAS 06530—0213, 07430+1115, and 
19477+2401) and the others are observed at higher S/N than in previous spectra. 
These have the same shape and central wavelength (20.1 //m) as found in the ISO 
spectra of the brighter PPNs. The 30 //m feature was seen in all seven objects. 
However, it is not resolved into two separate features (26 and 33 /xm) as was 
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claimed on the basis of ISO spectra, which presumably suffered from the noisy 
detector bands in this region. All showed the infrared aromatic bands (AIB) at 
11.3, 12.4, and 13.3 /im. Five of these also appear to have the C2H2 molecular 
band at 13.7 /xm, one in absorption and four in emission. This is extremely rare, 
with only one other evolved star, IRC+10216, in which C2H2 emission has been 
observed. Four also possessed a broad, unidentified emission feature at 15.8 /im 
that may possibly be related to the 21 /im feature. Model fits were made to 
the spectral energy distributions for these PPNs to determine properties of the 
detached circumstellar envelopes. The 21 /xm feature has been seen in all Galac- 
tic carbon-rich PPNs observed, and thus its carrier appears to be a common 
component of the outflow around these objects. 

Subject headings: circumstellar matter: — infrared: stars — ISM: planetary 
nebulae: general - radiative transfer — stars: AGB and post-AGB 



INTRODUCTION 



The study of infrared emission features in evolved stars provides an important avenue 
to study the chemistry of the circumstellar environment. First seen in the IRAS spectra of 
four sources (IKwok. Volk. fc Hrivnaklll989l ). a new, broad feature was discovered at ~21 /im 
(the so-called "21 /im feature"). Subsequently, a small number of additional sources with 
this feature have been identified based on spectra from IRAS, ISO, and with the CGS on 
the United Kingdom Infrared T elescope. The total number of such sources is about a dozen 
( iKwok. Volk. fc Hrivnaklll999al ). An ISO study of the strongest sources showed that they 
all had the sa me peak wavelength, 20.1 /im , and nearly the same shape with no evidence of 
substructure (IVolk. Kwok. fc Hrivnaklll999l ). 



The 21 /im sources all possess two characteristics in common. First, the 21 /im sources 



later confirmed in d etailed high-resolution abundance studies ( Van Winckel fc Reynier^ 



are all carbon rich. This was first deduced based on the presence of molecu. 
C3 in their visible spectra ( Hrivnak fc Kwok 1991a : Hrivnak 1995[ Bakker et al 



ar Co and 



19971) and 



2000 



Reddy et al.ll2002l ). The second common characteristic of the 21 /xm sources is that they 
are nearly all proto-planetary nebulae (PPNs), objects in transition between the asymptotic 
giant branch (AGB) and planetary nebula (PN) stages of the evolution of intermediate-mass 
stars. During this stage, the extensive AGB mass loss has ceased and the central star is 
surrounded by a detached, exp anding envelope of molecular gas and dust. This stage can 
last several thousands of years (jBlockerlll995bl ). during which the central star evolves from 
spectral type late-G to early-B, after which time it becomes hot enough to photoionize the 
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gas and pro duce a PN. (The 21 /xm feature also appe ars to be present but weak in a few 
young PNs; IHouv. Waters, fc TielenslboOll : IVolk Ibood l 



There has been an active search for the carrier of the 21 /im feature and many candi- 
dates have been suggested, including PAHs, fuUerenes, and nanodiamonds. Most but not all 
of the proposed carriers involve carbon. A good match to the feature is produced by TiC 
(Ivon Helden et al.ll200Cll ). but this requires severe formation constraints due to the low abun- 
dance of Ti, even in stars like these that have enhanced s-process elemen ts. SiC also provides 
a fit i n this region and without the severe abundance constraints (see ISpeck fc Hofmeister 
(120041 ) and references therein). However, most carriers suggested have associated emission 
features at other wavelengths that are absent or unconfirmed in these sources. 

A very broad, so-called "30 /im feature" ha s been know for several decades to exist in 
the spectra of evolved C-r ich AGB stars and PNs (IForrest. Houck. &: McCarthylll98ll ). It has 
been attributed to MgS (iGoebel fc Moseleylll985l ). which peaks in the range of 30—35 fim 



depend ing on the temperature of the grain. Subsequent stu dies have supported this identifi 



cation (iSzczerba et al.lll999l : iHony. Waters, fc Tielend 120021 ). It has been seen with ISO in a 



number of objec ts ranging from AGB stars to PPNs to PNs (IHony. Waters, fc TielensI 12002 



Volk et al.l |2002| ) on a dust continuum that peaks at monotonically longer wavelengths in 
the transition between these stages. On the basis of recent ISO spectra, it has been claimed 
that the feature has been resolved into two features ceii t ered at 26 and 33 /im in PPNs, PNs, 
and extreme- AGB stars jHrivnak. Volk. fc KwoM|200oI : IVolk et al.ll2002h . 



A family of emission features at 3.3, 6.2, 7.7, 11.3, 12.4, and 13.3 /im (historically known 
as unidentified infrared [UIR] bands) is commonly observed in the diffuse ISM. Although the 
specific carriers are not identified, it is quite certain that these arise from the stretching and 
bending modes of carbonaceous materials with an aromatic component, such as polycyclic 
aromatic hydrocarbons (PAHs). We will refer to these as aromatic infrared bands (AIBs). 
Since these AIBs are also observed in carbon-rich PNs and PPNs, the carriers of these bands 
are likely to be synthesized in the circumstellar environment of late-type stars. These are 
often accompanied in PPNs by features at 3.4 and 6.9 /im attributed to aliphatic groups 
attached to the aromatic hydrocarbons (jHrivnak. Geballe. fc Kwokl 120071 ) . The 11.3, 12.4, 
and 13.3 /im features lie within the wavelength range of the Spitzer spectrographs. 

In this study, we have obtained Spitzer spectra of seven PPNs to further study the 21 
/im, 30 fim, and AIB features. These targets were selected with the goal of completing the 
mid-infrared study of carbon-rich PPNs. In addition to new targets that met this criterion 
(IRAS 06530-0213, 07430+1115, 19477+2401), we included some carbon-rich PPNs known 
or suspected to possess the 21 /im feature based on low S/N ISO or ground-based spectra. 
Also included for comparison was one additional object that was well observed with ISO but 
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faint enough for Spitzer observations. 



OBSERVATION AND DATA REDUCTION 



The target objects were observed with t he Infrared Spectro graph (IRS; iHouck et al. 



20041 ) on the Spitzer Space Telescope {Spitzer^ IWerner et al.ll2004l ). Most were observed as 
part of program No. 20208, PI: B. Hrivnak. The targets were all relatively bright for Sptizer 
and consequently were all observed in the high resolution modes, Short-High (SH) and Long- 
High (LH). These modes consisted of ten overlapping orders, and ranged in wavelength from 
10—19.6 /im for SH to 19—37 fim for LH. The slit sizes for the two modes are 4.7" x 11.3" for 
SH and 11.1" x 22.3" for LH. Each observation consists of two nod positions, located at one- 
third and two-thirds of the way across the slit; at each nod position three scan cycles were 
carried out. The spectral resolution is R ~ 600. No off-source background measurements 
were made for these bright targets. An observing log is given in Table [TJ 

Two of the targets were initially observed using inaccurate coordinates and were re- 
observed during Director's Discretionary time. Three of the objects were not observed in this 
program in the SH mode because they were on the GTO program no. 93, PI: D. Cruikshank. 
Data for these were obtained from the Spitzer archive and have been incorporated into this 
study. Their observations were made in a similar manner to those described above. Details 
of these additional observations have also been included in Table [H 

Unfortunately, the coordinates of some of the targets were not know accurately and 
thus the stars were not always centered well on the slits. Accurate positions are important 
since the slits are relatively small compared to those used in previous mid-infrared spectral 
observations. In some cases, this meant that the spectrum at one or both of the two nod 
positions was weak or noisy or could not be extracted properly. Peak up was performed with 
nearby stars using the PCRS. To avoid such problems in follow-up studies, accurate positions 
are list in Table [2] for each of the targets as derived from accurate (<1") 2MASS positions. 
For the fainter visible sources, source identifications are based on finding charts determined 
from our ground-based mi d-infrared observation s. Fi nding charts for IRAS 194 77+241 l|l| and 
22574 -f6609 are shown by lHrivnak et al.l (119851 ) and lHrivnak &: Kwokl (jl991bl ). respectively. 



The data were initially processed using the data reduction pipehne at the Spitzer Science 



^ Based on comparis on with the 2MA SS catalog, it appears that the JHKL' magnitudes that we pubUshed 
for IRAS 19477+2401 |Su et alllioOll ) were for the nearby source located 8" west and 5" south, 2MASS 
19495429+2408481. The 2MASS magnitudes for IRAS 19447+2401 are J=12.61, H=10.75, and K=9.61. 
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Center (SSC): version 13.2.0 for the earlier data and ver sion 16.1.0 for the two that were 



re-observed. We then used the software program SMART (iHigdon et aLll2004l ) to reduce the 
data. We started with the Basic Cahbration Data (BCD) files. These were first cleaned for 
bad pixels using IRSclean. Bad pixels included the so-called rogue pixels, which are identified 
by observing campaign at the SSC, and other bad data which are flagged in the pipe-line. The 
spectra were then extracted and defringed. The resulting spectra were examined manually, 
by order, within each of the individual scan cycles. Bad data were removed from the edges 
of the orders. Data spikes not present in all the individual scan cycles or in both of the 
overlap regions in adjacent orders or in both of the nod positions were treated as artifacts, 
as were some additional large spikes, and these were removed. An average spectrum was 
composed in each nod position. In some cases we made small shifts in flux (2—10%) to bring 
the overlapping orders into better agreement. The spectra from the two nod positions were 
then compared and an average spectrum formed if the two spectra agreed. In some cases, as 
mentioned above, the spectrum from one of the nod positions was near the edge of the slit 
and the data were of inferior quality; these spectra were not included. The spectral region 
14.5—14.6 /um had lots of spikes in the data and was not a region with overlapping orders; 
this resulted in particularly noisy spectra in this region. The LH spectra beyond 29 jjm were 
more noisy and those beyond 36 fim were completely deleted. 

The SH and LH spectra for each object were then combined, scaling where necessary 
based on a comparison of the region of spectral overlap from 19.0 to 19.6 ^m. In all cases the 
SH was fainter or similar to the LH; this may be due to the larger slit size of the LH mode, 
better accommodating small pointing errors. Therefore we scaled the SH when necessary 
to match the LH in the overlap region to produce a final Spitzer spectrum for reach object 
covering the region 10 to 36 fim. 

We note for comparison purposes that previous mid-infrared spectra have been obtained 
of some of the target objects. Initial observations were at low resolution from space with LRS 
(7.7—23 /im, R~25) on IRAS, or from the air using the Kuiper Airborne Observatory (KAO; 
16—48 ^m, R~33), or from the ground using CGS3 (7.5—13 ^m, R~50; 16—24 /im, R~72) 
on the 3.8-m United Kingdom Infrared Telescope (UKIRT). Higher resolution observations 
were made more recently with SWSOl (R~250) or SWS06 (R~600-2000) on ISO. Some 
of these will be referenced in our following discussion of the new higher-resolution, higher 
signal-to-noise ratio Spitzer spectra. 



3. THE SPITZER IRS SPECTRA 



IRAS 23304+6147: The SH and LH spectra in the two nod positions were in good 
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agreement and were averaged. In the LH, the longest wavelength order (LHll: 34—36 fim) 
was scaled upward by 5%. The complete spectrum is shown in Figure [H The spectrum 
shows strong broad 21 fim and 30 fim emission features. Also seen are medium-width 
emission features at 11.4, 12.3, 13.2 /xm, a medium-width feature at 15.9 fim and a weak 
feature at 22.3 /zm. The short wavelength SH region is shown separately in Figure [2], where 
the features in that wavelength region can be seen more clearly. 

IRAS 23304+6147 is a well-documented 21 ii m source, being one of the ori ginal four 
IRAS sources in which this f eature was detected (IKwok. Volk. fc HrivnakI Il989l ). The 30 
fim feature was discovered bv lOmont et al.l (119951) using the KAO. Good observat ions were 
subsequently carried out with ISO (IVolk. Kwok. &: HrivnakI Il999l : I Volk et al.ll2002l ). and we 
observed it again partly for comparison purposes. Shown in Figure [3] are the ISO SWSOl 
(3—45 /im) and PHOT-S (3—8 fim) and the KAO (16—48 fim) spectra along with the Spitzer 
spectrum. It can be seen that the spectra are very similar in intensity from 10—19.5 /im and 
similar in shape from 10—27.5 /im. Beyond this, out to 36 /im, the ISO spectrum is about 
10% higher. The KAO spectrum is roughly similar in flux to the Spitzer spectrum beyond 27 
/im, but neither it nor the ISO spectrum shows the flattening seen in the Spitzer spectrum 
beyond 33 /im. The new Spitzer spectrum is noticeably fiat between 30 and 35 fim. 

The Spitzer spectrum was convolved with the IRAS bands to calculate synthetic IRAS 
photometry. This resulted in a 25 /im flux of 59.9 Jy, in excellent agreement with the IRAS 
25 /xm flux of 59.1 ± 4% (53.7 Jy color-corrected). The color-corrected IRAS 12 and 25 /im 
flux measurements have also been included for comparison in Figure [31 

IRAS 05113+1347: The spectra in each of the two nod positions in SH and LH agreed 
with each other and were averaged. To combine SH and LH smoothly required that we 
increase the SH spectra by 7%. The spectrum shows a moderately strong 21 /xm feature and 
a strong 30 /im feature. At the shorter wavelengths are seen emission features at 11.4 and 
12.1 /im and perhaps weak features at 13.3 and 15.8 ^m, along with a weak feature at 22.3 
/im. 

The 21 //m feature was flrst detected in this object i n a low S/N CGS3 spectrum from 
UKIRT following the suggestion of it in the IRAS spectra (IKwok. Hrivnak. &: Geballdll995l ). 
but the source was not observed with ISO. The Spitzer synthetic IRAS photometry resulted 
in a 25 /zm flux of 14.2 Jy, 7% lower than the IRAS 25 /im flux of 15.3 ± 6% (12.6 Jy 
color-corrected) . 

IRAS 05341+0852: The SH and LH spectra in the two nod positions agreed reasonably 
well with each other and were averaged. The flux levels in most of the SH orders were scaled 
by a few percent to agree with the adjacent orders, but the data for the flrst three orders 
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(10—12 fim) agreed well without any shifts. The SH and LH spectra were then combined 
without any shift in flux. The spectrum shows a strong 11.3 /im feature, a moderately strong 
21 fim feature, and a strong 30 feature. 

The 21 nm feature in this so urce was also first detected in CGS3 spectra from UKIRT 



( iKwok. Hrivnak. fc Geballdll995l ): the source was not observed with ISO. Present is a fea- 
ture at 12.4 fim, and a broad but weak feature at ~15.8 fim. Also seen is a very narrow 
emission feature at 13.7 /im that agrees with the position of the C2H2 molecular band seen 
in absorption in carbon stars. The simulated IRAS photometry value was found to be 9.4 
Jy at 25 /im, compared to the observed IRAS value of 9.9 ± 7% Jy (8.7 Jy color-corrected). 

IRAS 06530—0213: The spectrum was good at only one nod position in each of SH 
and LH. The longest wavelength order in LH was scaled up by 8% to match the rest of 
the LH spectrum. The SH spectrum was scaled up by 24% to merge the two spectra; 
this is presumably to correct for flux missed due to poor centering of the SH spectrum. 
(Note that if poor centering is indeed the cause of the flux loss, the proper correction is 
wavelength dependent and not simply a constant scale factor. However, we do not have 
enough information to make such a correction; thus this likely results in an uncertainty of 
several percent in the shape of the SH spectrum.) The 21 /im feature is strong; the 30 /im 
feature is not as relatively strong as in IRAS 23304+6147 and it shows an unusual inflection 
around 33 /im. Also seen are AIB emission features at 11.4 (strong), 12.3, and 13.2 /im, 
a strong, medium-width feature at 15.8 /im, and a weak feature at 22.3 /im; the apparent 
feature at 14.6 /im is likely an artifact. There appears to be a weak, very narrow emission 
feature at 13.7 /zm, the position of the molecular C2H2 band. 

This object had no previous mid-infrared spectrum, so this is the discovery spectrum of 
the 21 and 30 /im emission features in this object. The Spitzer?lxD^ is noticeably lower than 
the IRAS ?l\xyi, with a synthetic 25 /im flux of 21 Jy as compared with the IRAS?l\x:^ of 27.4 
± 4% Jy (25.7 Jy color-corrected). This large difference in flux is attributed to flux loss in 
the Spitzer spectrum due to the non-optimal centering of the object in the slit. 

IRAS 07430+1115: In both the SH and LH spectra the two nod positions are in good 
agreement and were averaged. The SH and LH spectra were combined with no shift in flux, 
but the very shortest part of the LH spectra was removed due to lack of agreement with the 
SH in part of the region of overlap. Seen is a 21 /im feature and a strong 30 /im feature. The 
11.3 /im feature is seen a broad plateau extending out to 12 /xm. The Spitzer flux is close to 
the IRAS flux, with a synthetic 25 /im flux of 27.9 Jy, close to the IRAS flux of 29.9 ± 5% 
Jy (24.7 Jy color-corrected). 



The earlier UKIRT CGS3 spectrum did not reveal the 21 fim feature 



Hrivnak Sz Kwok 
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( 119991 ). but since the CGS3 spectrum only extended to 24 fim and at lower S/N, it is not 
surprising that this weak feature was not detected. This target was not observed with ISO. 
Thus this is the discovery spectrum of both the 21 and 30 //m features. Based on the earher 
CGS3 spectrum, it had appeared that IRAS 07430+1115 was an exception to the rule that 
all C-rich PPNs are 21 fim sources. This exception is now removed. 

IRAS I9477+24OI: Only one of the LH nod positions contained a good spectrum. Nei- 
ther of the SH nods had a useable spectrum; it was on the edge of the slit and appeared very 
faint. These problems were traced to poor positioning of the object in our spectrum and 
mis-identification of the source in the SH archival observation. The originally published coor- 
dinat es were not of sufficient accuracy, but a correct finding chart was shown (IHrivnak et al. 
19851 ). In the one LH spectrum that we used, there was still a need to scale some of the 
individual orders by a few percent before combining them. 

The resulting spectrum is shown in Figure 1. It shows 21 /xm feature and a strong 30 fim 
feature. In addition, there is a weak feature emission at 22.3 /im and possible weak emission 
features at 28.6, 29.8, 32.6, and perhaps even 34.9 fim. At 31.4 yum there is a narrow dip 
followed by a narrow peak; since it appeared in two separate orders, we did not remove it 
from the spectrum. However, it would be quite unusual to see a P Cygni profile at this 
spectral resolution. The 30 //m feature shows a slight infiection in the spectrum around 33 
/xm. The synthetic IRAS 25 /im fiux is 50.5 Jy, 8% less than the IRAS value of 54.9 ±4% 
(51.2 Jy color-corrected). 

IRAS 19477+2401 was observed with CGS3 but the 2 1 /im feature was not seen in 



the noisy 20 fim spectrum (IKwok. Hrivnak. fc Geballe 



1995r). and the ISO spectrum shows 



fiux problems also attributed to an inaccurate position (IHrivnak. Volk. fc Kwokll2000l ). This 
Spitzer spectrum is the discovery spectrum for the 21 /im and the 30 fiin emission features 
in this source. 

IRAS 22574+6609: In both SH and LH, the spectra in the two nod positions agree well 
and were averaged. In the LH, we shifted the longest order (LHll) up by 6% to match with 
the previous order, and the SH was scaled up by 2% to match with the LH. The spectrum 
shows a moderately strong emission features at 21 /im and 30 /im. A very strong emission 
feature is seen at 11.3 /xm, along with a medium-width feature peaking at 12.6 /tm and 
another one at 17 fim (about 0.5 /im wide). The is evidence of a weak absorption feature at 
13.7 /im that may be the molecular band due to C2H2. 



The 21 /im feature was first detected in this object in the IRAS LKS spectrum (IHrivnak fc Kwok 



991br). followed by observati on of the 21 and 30 um and AIB features with ISO SWSOl 



(IHrivnak. Volk. fc Kwokll2000[ ) and SWS06 spectra (IVolk. Kwok. fc Hrivnaklll999l : IVolk et al 



- 9 - 



20021 ). These previous spectra were rather noisy and the new Spitzer spectra show the 21 



and 11.3 fim features at much higher S/N. In Figure [3] are compared the Spitzer and ISO 
spectra. They have the same general shape out to 35 //m, but the Spitzer spectrum is ~20% 
higher shortward of 23 /im. The spectral features differ in the region 12 to 14 /im; the Spitzer 
spectrum shows a feature at 12.6 fim while the noisier ISO spectrum shows features at 12.0 
and 13.3 /um. The Spitzer synthetic 25 fim flux is again less than the IRAS flux, 27.3 Jy 
compared with 29.5 ±4% Jy (29.8 Jy color-corrected). 

We note that in most of the sources (IRAS 06530-0213, 07430+1115, 19477+2401, 
22574+6609, 23304+6147) the spectra are flat or even show an upturn in flux beyond ~33 
//m. In the previous spectra of IRAS 23304+6147 and 22574+6609 this was not seen. As 
mentioned in the discussion of the individual sources, the longest wavelength order, LHll 
(34—36 /um), was scaled upward by 5—8 % of some cases to match the rest of the LH 
spectrum. Although the data are of high S/N, we maintain some skepticism about the 
accuracy of the shape of the spectra in the 33—36 fim region. Perhaps slit losses make some 
contribution the the suspected problem, since some of the sources are know to be extended 
in the N band. 



4. DISCUSSION OF THE SPECTRAL FEATURES 
4.1. 21 /im Emission Feature 



Volk. Kwok. fc HrivnakI (119991 ) used ISO spectra to examine the 21 fim feature in eight 
PPNs, including IRAS 23304+6147 and 22574+6609. They found the feature to have a peak 
wavelength of 20.1 ± 0.1 /im with a rather similar shape in all eight sources. However, they 
found a large range in the relative strength of the feature (peak-to-continuum ratio). 

We can similarly compare the feature profile in these new, higher S/N Spitzer spectra. 
This was done for the six objects with spectra in both SH and LH modes (for continuum 
fitting). To determine the continuum, we used the flux at the shorter and longer wavelength 
sides and fit it with a spline function that gave a good smooth fit to this continuum; these 
ranged in order from 5 to 8. For three of the sources, IRAS 06530—0213, 22574+6609, and 
23304+6147, we treated as the continuum the flux at 13—18 and at 24—33 /im (which 
also includes the part of the 30 fim feature), while for the other three, for which there was a 
more steeply rising continuum on the long-wavelength side, we used a smaller range about the 
feature. This function was then divided into the spectrum to produce the feature profile. The 
relative strength of the feature varied quite a bit, from 0.21 to 1.46. An approximately similar 



range was determined in previous studies of the brighter objects ( jVolk. Kwok. &: Hrivnak 
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19991 : iHrivnak. Volk. fc Kwoki l200d : IVolk et al.l |2002| ) ; however these present spectra are 
superior due to their higher S/N. These peak-to-continuum ratios are hsted in Table [3l 
along with the feature width and peak wavelength for each source. The feature profiles are 
similar; this is shown graphically in Figure HI where the profiles have been normalized to 
1.000 and the continuum removed. The long wavelength side of the feature differs slightly 
among the sources, especially for IRAS 05113-1-1347. However, this depends sensitively on 
the 30 /im feature, and the assumption that there is a region of continuum between the two 
features and thus that the 30 /xm feature does not contribute flux between 21 and 24 /im. 
The short wavelength side of the 21 /im feature is less susceptible to contamination by strong 
features and appears to be relatively consistent among the sources. The peak wavelength of 
the feature is 20.1±0.1 /im, the same as found for the sources in the earlier ISO study. The 
feature width we find is greater than in the earlier study, but this is likely due to the fact 
that we had more continuum to work with. IRAS 23304-1-6147 was included in both studies, 
and its 21 /im feature is wider by 20% in the Spitzer spectrum. 

In all of the 21 /im sources studied in detail, this feature appears to have a relatively 
similar central wavelength and profile shape in each, but different strengths. Uncertainties 
arise in defining the feature since it sits upon a continuum and adjacent to a very broad 30 
/im feature, both of which are rising steeply toward longer wavelengths. 

All carbon-rich PPNs that have been observed possess the 21 /tm emission features 
jHrivnak et ahlboosh . The sole exception to this statement may be IRAS 01005-f7910 (OBe). 
The low S/N ISO spectrum of this s ource may or may not possess the 21 /im feature; it does 
possess the 30 /im and AIB features (IHrivnak. Volk. fc Kwokll2000l ). Since this is the hottest 
know carbon-rich PPN and since the feature appears to be absent or weak in PNs, it may be 
that IRAS 01005-1-7910 is displaying the expected weakening of this feature with evolution 
toward a hotter central star. The next hottest carbon-rich PPN star is IRAS 16594—4656 



(B7e), and it clearly shows the feature (IGarcia-Lario et al.l 119991 ). 



The carrier of this feature is still not identified. There have been many suggestions, 
inclu ding PAHs, hydrogenat ed fuUerenes, nano diamonds, and more recently TiC naii oclus- 



ters (Ivon Helden et al.ll2000l ). doped SiC (ISpeck fc Hofmeistei 



the interaction of Ti atoms with fuUerenes (iKimura et al. 



2004J : I Jiang et al.ll2005f ). and 



20051). (For a review of these and 



other potential carriers and additional references, see IPosch et al. 



2004J ) . Some of these pre- 



dict features at other wavelengths, but none of them have been confirmed. SiC does indeed 
produce an obs erved feature at 1 1.3 /im, but it does not have the correct flux ratio with the 
21 /im feature (I Jiang et al.ll2005l ). We note below a possibly correlated feature at 15.8 /im 
in our observed spectra. 
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4.2. 30 fim Emission Feature: Two Components? 



We have previousl y used ISO spectra to study the 30 fim feature in PPNs ( iHrivnak. Volk. fc Kwok 



2000l : IVolk et al.ll2002l ). In these ISO spectra it appeared that this very broad feature was 



resolved into two separate features at 26 and 33 /im, seen with varying strengths in PPNs. 
These new Spitzer spectra do not support this claim. This can be seen by the comparison 
of the Spitzer spectrum with the ISO spectrum for IRAS 23304+6147 as shown in Figure 
[3^. The apparent downward inflection in the ISO spectrum at 27.5—29.5 /im is evidently 
due completely to the unreliable flux of ISO band 3E (27.5—29.5 /im). Also, the data in 
ISO band 4 (29.5—45 /im) are known to be very noisy due to the many cosmic ray hits 
and the flux calibration is less reliable. A similar effect is seen in the comparison of the 
Spitzer spectrum of IRAS 22574+6609 with its ISO spectrum as shown in Figure [3]d. The 
downward inflection in the ISO spectrum between 26 and 33 /xm is not seen in the Spitzer 
spectrum of IRAS 22574+6609, nor is it seen in any of our other Spitzer spectra of sources 
with the 30 /im feature. Low-resolution Spitzer spectra of evolved c arbon stars in the LMC 
also do not show an indication of the 30 fim feature being resolved (iBuchanan et al.ll2006l ). 
although this feature does appear to be significantly weake r in these LMC sources than in 
similar objects in the Milky Way Galaxy or in these PPNs. iHony. Waters, fc Tielend (120021 ) 
used ISO spectra to study the 30 fim feature in a large sample (63) of carbon-rich evolved 
stars — AGB, post-AGB, and PN. They treated it as a single feature and found that a good 
fit to the 30 /zm profile could be obtained using MgS grains with a temperature differing 
from the continuum temperature. However, in ~40% of these a residual excess existed at 26 
/im. Whether this excess is also an artifact due to the above problems in the ISO spectra 
might be further investig ated using Spitzer spectra of these ISO targets. In a recent paper. 



Zhukovska &: Gail I (120081 ) present a strong case for the feature to be due to MgS as a coating 



on SiC grains. They go on to show how the feature shape and peak would change with the 
relative thickness of the coating, with the main peak at 25—28 /xm and predicting a sub-peak 
at 33—38 /im when the coating is not too thick. 



4.3. AIB Emission Features 

A strong 11.3 /im feature, along with several corresponding features in the 12—14 /im 
region and an underlying plateau. The known features at 11.3, 12.1, 12.4, and 13.3 /im are 
due to out-of-plane bending modes of aromatic compounds with peripheral H atoms attached 
to the aromatic rings. Although the carriers of these features are frequently attributed to 
PAH molecules, the strong associated underlying broad continuum suggests that the carrier 
of the continuum emission cannot arise from small, gas-phase molecules, and has to be 



- 12 - 



emitted by solid-state grains or very large molecules with thousands of C atoms. 



To investigate the shapes of these features, we show in Figure O the 10—17 //m spectral 
region with the continuum removed (by division of a fourth-order cubic spline fitted from 
10.0—10.5 /im to 14—18 /im). The strong 11.3 /im feature has a relatively consistent peak 
wavelength of 11.33 ± 0.03 /im (ranges from 11.27 to 11.36 /im in these six PPNs), Shown for 
comparison is the ISO spectrum of IRAS 21282+5050, a young PN with a [WC] central star. 
Its spectrum in this region has a general similarity to that of the PPNs. However, the 11.3 
/im feature is narrower and peaks at a shorter wavelength, 11.21 /tm. The peak of the 11.3 
/im feature at shorter wavelengths is what is also seen in most other sour ces. From a s ample 



of 15 objects, consisting of HII regions, young stellar object s, and PNs. iHonv et al 



found that they all had a profile peak very near to 11.23 //m. Ivan Diedenhoven et al. 



( 2001 ) 

JiooJ 



defined a classification scheme for the 11.3 //m feature that depended upon the shape and 
peak wavelength of the feature, but all three of their clas s es (A, AB, B) had peak wavelengths 
ranging from 11.20 to 11.25 /im. However, ISloan et al.l (120071 ) found in a sample of objects 
with cooler central stars that the central wavelengths are shifted to the red relative to the 
above peaks found for the class A and B so urces. They include class C AIB features, based 



Peeters et al.l (120021 ). and these peak at wavelengths similar to what 



on the c lassification of 

we find. Sloan et al. J2007h attribute the difference in the class C AIB spectra as due to the 



fact that "the carbonaceous material has not been subject to a strong ultraviolet radiation 
field, thus allowing the relatively fragile aliphatic materials to survive." Our targets all 
have cool central stars and aliphatic features (see below), and thus are consistent with this 
suggestion. 

Around 12 /im, a broader emission feature is seen in all of these spectra, with a 
peak wavelength ranging from 12.0 to 12.6 /im. In IRAS 06530-0213 and 23304+6147 
and perhaps 05113+1347, a weaker emission feature is seen at 13.2—13.3 /im. The fea- 
tures at 12.1, 12.4, and 13.3 /im have been identified as the duo, trio, and quatro out-of- 
plane bending modes of C—H bonds and are associated with the 11.3 //m solo-CH feature 
( iHudgins &: Allamandolalll999l ). In previous studies of the infrared spectral characteristics 
of PPNs, we fou nd that unlike PNs and HII reg ions, PPNs show aliphatic emission fea- 
tures at 3.4 /im ( iHrivnak. Geballe. &: Kwokl 120071 ) and 6.9, in addition to the features at 
12.1, 12.4 and 13.3 /im which are uncommon in PNs or in other galactic nebulae. This 
spectral difference is suggested to be the result of chemical evolution, with a progressive 
formation of larg er and larger clusters of arora atic rings and the removal of peripheral 
aliphatic groups (IKwok. Volk. fc HrivnakI Il999bl ). Absent in these PPNs is the relatively 
strong, asymmetric 12.7 /im feature found in many objects with AIB features, including 
IRAS 21282+5050 JHouv et al.l I2OO1I : Isioan et all I2OO5I . I2OO7I 1 In this sample, it is found 
only in IRAS 06530-0213. 
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The existence of broad plateau emission features at 8 and 12 ix m are also a common fea- 
ture a mong; the PPNs but not PNs and objects with hot central stars (IHony. Waters, fc Tielens 



200ll ). This emission plateau is attributed to a collection of out-of-plane ben ding modes of a 



miscellaneous mixture of aliphatic side groups attached to the aromatic rings (iKwok. Volk. fc Bernath 



20011 ). From Figure [5l we can see that the 12 /im plateau feature is present among most of 



the PPNs observed here, with IRAS 05113+1347 perhaps being an exception. The width 
of this feature is estimated to be about 2 fim. An exception here is IRAS 07430+1115, 
which appears to have a broader width extending to 10 /um, although it is hard to set the 
continuum at the short end of this spectrum. 



4.4. 15.8 /im Emission Feature 

A relatively strong and medium-wide (~1.3 /im) feature is seen at 15.8 /xm in two of 
the sources, IRAS 23304+6147 and 06530—0213. The feature a l so ap pears to be present 



but weak in IRAS 05113+1347 and 05341+0852. iMoutou et al.l (|l999r ) point out features 
at 15.8 and 16.4 fim in their spectrum of NGC 7023, noting that their widths of ~0.15 fim 
are similar to that of the 11.3 fim feature. In the more e x tensiv e discussion of sources with 



features between 15 and 20 fim by Van Kerckhoven et all (2000), this feature can be seen in 



two or three sources but is not particularly distinct and is usually weaker than the apparently 
associated feature at 16.4 /fm, which is not present in our spectra. These and the complex 
of other features which are thought to make up a 15—20 /xm emission plateau in carbon-rich 



sources are attributed to C-C-C bending modes in PAHs. ISellgren. Uchida. fc Werner! (120071 ) 
point out additional emission features in NGC 7023 at 17.4 and 17.8 /tm, along with a broad 
plateau feature at 17 /im and a weaker features at 18.9 /im. Our C-rich PPNs show 3.3, 6.2, 
11.3 /im features similar to those seen in NGC 7023 but none of these other, new features. 
The 15.8 fim features seen in their spectra are not nearly as broad or strong as the 15.8 fim 
features seen in several of our new PPN spectra, and thus they are not the same. 

The width of this feature suggests that it is made up of a larger molecule than that 
producing the AIB features and perhaps is due to a solid-state feature. To quantify the 
properties of this feature, we determined it peak wavelength and the relative strength of the 
feature (peak-to-continuum ratio). For the weaker sources, we determined only an upper 
limit to the relative strength. The values are as follows: IRAS 06530—0213 (15.75 /im, 0.22), 
23304+6147 (15.79 /im, 0.12), 05341+0852 (15.77 /im, 0.06), 05115+1347 (..., <~0.03), 
22574+6609 (..., <-0.03), and 07430+1115 (..., <0.02). In the source with the strongest 
feature, IRAS 06530—0213, it has a width of 1.3 fim (15.1 to 16.4 fim). The shape in IRAS 
23304+6147 is consistent with this and the shape in the weaker sources is also consistent. 
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but the feature in these is of much lower S/N. 

As can be seen by comparison with Table [31 the two sources with the much stronger 15.8 
fim feature are also the sources in our study with the much stronger 21 /im feature. This 
raises the possibility of a possible relationship between these two features. An examination 
of the lower S/N ISO spectra of other 21 /im sources shows that the 15.8 /im feature is 
likely present in IRAS 22272+5435, 19500-1709, 16594-4656, 23304+6147 (as confirmed in 
this Spitzer spectrum), and in the low S/N spectrum of IRAS 07134+1005, the sourc e with 
the relatively brightest 21 /im feature JVolk et allbooi iHrivnak. Volk. fc Kwoklboooh . This 
possible relationship is based on only a few objects and so is only tentative, but it would 
certainly be worthwhile to investigate it for more of the 21 /im sources. The presence of a 
feature correlated with the strength of the 21 /im feature would provide an important clue 
to help identify the carrier and could at minimum help to constrain the viable suggestions. 



4.5. 13.7 //m C2H2 Feature 



We have previously mentioned a narrow (unresolved) feature at 13.7 /im in three of these 
PPNs, in absorption in one (IRAS 22574+6609) and in emission in two (IRAS 05341+0852 
and 06530—0213). This is at the wavelength of the m olecular C^^H^ (acetylene) vb bandhead 



feature commonly seen in ab sorption in carbon stars flWillemslll988l : ICernicharo et al.lll999 



Volk. Xiong. fc Kwok 1120001 1. This feature may also be present in IRAS 23304+6147 in both 
absorption and emission (P Cygni profile; see Figure [5]) and in IRAS 07430+1115 in emission. 
This C2H2 f eature has previously b een reported in absorption in only three galactic PPN, 



AFG L 618 flCernicharo et al.ll200lh . AFGL 2688, and IRAS 13416-6243 flKraemer et al. 



2006h: it has recently been also seen in absorption in the LMC PPN/P N object SMP LMC 1 1 



(Bernard-Salas et al 



20061 ) and the SMC PPN object MSX SMC 029 flKraemer et al.ll2006h . 



Kraemer et al.l (l2006l ) and ISloan et al.l (l2007l ) point out that when the C2H2 absorption 
feature appears in so urces with PAH era ission, the sources are those classified as the rare 
class C PAH sources ( iPeeters et al.ll2002l ). Our sources with C2H2 would all be classified as 
class C PAH sources based on their broad 8 /im features and secondarily based on the longer 
wavelength of their 11.3 /im feature, and thus this claimed correlation is supported by our 
new observations. 

To our knowledge, C2H2 has pr eviously been seen in emi ssion in only one evolved stellar 
object, the AGS star IRC+10216 JXeadv fc Ridgwavl E993[ ) . Thus this is the first report 
of the molecule in emission in a post-AGB object. One might consider the emission to be 
due to collisional excitation. However, co llisionally-excited H2 em ission is not seen in any of 
these three PPNs at 2.12 /im (1-0 S(l)) (IKelly &: Hrivnakll2005l ). nor is H2 emission seen in 
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any objects in this study at 12.27 /im (0-0 S(2)) or 17.03 //m (0-0 S(l)). (Only for IRAS 
22574+6609, with a 13.7 /im absorption feature, is a (weak) 2.12 /im H2 emission feature 
seen.) We have recently confirmed the presence of C2H2 z^5 lines having P-Cygni profiles in 
three of these objects using high resolution spectroscopy, as will be discussed in more detail 
elsewhere. This feature has the potential to serve as a diagnostic of the physical conditions 
in the circumstellar envelope. 



4.6. Summary of Spectral Features 

We summarize the presence of these spectral features in our targets in Table HI Also 
included are the carbon-to-oxygen (C/0) abundance ratio and observations of AIB features 
at shorter wavelengths as found in other studies. It can be seen that all of the sources in 
this sample with the 21 /im feature are carbon-rich, possess the 30 /im emission feature, and 
also possess AIB emission features in the 3—14 /zm range (except IRAS 19477+2401, which 
has not been observed in high S/N spectra in this range). These spectral features appear 
to accompany the presence of the 21 /im feature. All of the sources in this study that have 
visible spectral types also possess C2 and in most cases C3 in absorption. In addition to the 
in-some-cases-strong emission feature at 15.8 /im, there is an unidentified weak feature seen 
at 22.3 /im in four of these sources. 



MODELING OF THE SEDS 



To learn more about the physical conditions under which these features arise (T, p), 
and also to determine the mass loss rates (M) and expansion time scales (tdyn), we carried 
out modeling of the circumstellar envelopes. Model fits to the spectral energy distributions 
(SEDs) were carried out for the six sources with complete 10—36 /im spectra. The contin- 
uum was assumed to be due to amorphous carbon (AC) grains, with optical constants by 
Rouleau fc Martini (119911 . type AC2) and an assumed size of 0.1 /im. The radiation transfer 
code by iLeungj (Il976l ). DUSTCD, was used to produce one-dimensional models. The dust 
envelope is assumed to have an inner radius r^^, a dust density profile of 1/r", and to ex- 
tend to 0.75 pc (somewhat arbitrarily set at the distance reached in 5 x 10^ yrs with an 
expansion rate of 15 km s~^). The temperature of the central star (T*) was initially set at 
the value derived in published atmospheric model analyses and then adj usted slightly (<250 



199 



i) with 



K) to fit the SED at shorter wavelengths. A Kurucz model atmosphere (iKuruca . 
log g = 1.0 was used for the central star, and this was the assumed energy source for the 
nebula. The interstellar extinction Ay was treated as an adjustable parameter. The other 
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parameters that were initially adjusted were the temperature (Td) and inner radius of the 
cool dust shell, the optical depth (tq), and the exponent of the density law {a). These were 
adjusted by trial-and-error fits of the calculated to the observed SEDs until a "best fit" was 
determined. We began by fitting the mid-infrared continuum at 18 /im and 100 and/or 60 
lim; we set to the highest value that did not go above the continuum and adjusted rin 
and a. We then adjusted the optical depth in the near-infrared, along with Ay and T*, to fit 
the near-infrared and also the visible as well as possible. These parameters are not strongly 
coupled, with each primarily effecting a dist i nct as pect of the SED. More details of the basic 



model are given by iHrivnak. Volk. fc Kwoki (120001 ) 



After initial attempts with this basic model to fit the SEDs of these PPNs, the model 
was adjusted in two ways. First, in some cases it was found that a single dust shell with the 
high density needed to fit the mid-infrared SED and a large radial extent produced flux that 
was too high to fit the long wavelength 60 and 100 /im IRAS finxes, even if we assumed a very 
steep density law (a>4). Thus we added to the model a discreet thin, dense, cool (c) dust 
shell of inner and outer radius Tin{c) and Tout{c), respectively. Second, it was found that for 
most of the sources, the model produced flux that was too low to fit the observations in the 
M and L bands. To improve this fit, a component of hot dust was added where needed. This 
hot (h) component was of low density and mass, and served to bring the model into better 
agreement with the observed flux in the 3—5 /xm region. The hot dust component begins 
close to the star, with a small value for rj„(h). The values for Trf(h) and rin(h) were adjusted 
to fit the M and L observations. The density laws for both the cool and hot components 
are assumed to have the same exponent. The resulting density ratio of the cool to the hot 
component (p(c)/p(h)) ranges among the models from 7700 to 15,000 at a radius of rin(c), 
and the hot dust contributes only 1.5—11 % of the optical depth at 11.2 /im. 

The AIB plateau features were not included in the modeling. They are likely produced 
by some combination of fiourescence and non-equilibrium heating, which the dust radiative 
transfer code does not have the capability to include. However, they do not carry a large 
part of the total energy emitted in the mid-infrared and presumably therefore would not 
have much of an effect on the optical depth in the dust shell or the derived parameters. 
IRAS 22574-1-6609 is a possible exception to this generalization since it has an unusually 
high optical depth dust shell. 

Empirical fits to the 21 and 30 /um features were added to the continuum. The 21 //m 
feature template is the same one we used previously to fit ISO spectra, derived from the 
ISO high-resolution SWS06 spectra of carbon-rich PPNs. The 30 /im feature template is 
derived from the fitting of the IRS spectrum for IRAS 05113-1-1347, since this source has a 
strong and relatively well-defined 30 /tm feature. The long wavelength end of the feature 
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was defined by a linear extrapolation of the observed feature shape from 32 to 36 /im, which 
goes to zero at 43 u m; that is reasonably consistent with the ISO l ong wavelength limit of 
the feature in PPN JVolk et al.lbood : IHouv. Waters, fc Tielei3 12002[ ) . (The extent to which 
this template fits or does not fit for the other objects is either telling us something about 
problems in the IRS data at long wavelengths or is indicating that there are variations in 
the feature shape.) We did not attempt to empirically fit the AIB plateau features. 

These model results are listed in Tableland the fits are shown in Figured In calculating 
the mass loss rate, a gas-to-dust density ratio of 330 was assumed. Note that while the total 
mass of the enhanced shell is well determined, the mass loss rate is less so, since a similar 
model fit could probably be obtained using a narrower or wider shell of the same mass. 
Below we discuss briefly the model for each object separately. 

IRAS 05113+1347. — The fit to the SED of this object required a discrete cool shell. 
Initial attempts to fit it with an extended dust shell required a very high density exponent of 
a = 5, and even then the fit was not as good as with a discrete shell. The hot dust component 
is needed to fit the M-band flux. This hot component is of low mass, and contributes only 
3% of the optical depth at 11.2 /im. The model provides a good fit to the SED, except 
for the emission plateau from 8 to 14 /zm and the 15.8 fim emission feature which were not 
included in the modelling. Th e visible image sh ows only the star, suggests that we may be 
looking pole-on to the nebula (lUeta et al.ll2000l ). 



IRAS 0534-1+0852. — The SED of this object was fitted by an extended, detached cool 
dust shell with a high density exponent of a = 3.5. To fit the L and M flux, a hot dust 
component is also needed. This hot component contributes 11% of the optical depth at 11.2 
/xm but less than 0.03% of the total mass in the envelope. The model fit to the SED is good, 
except for the emission plateau from 8 to 14 /xm and some excess emission observed from 
30 to 36 /im. The hot dust compon ent leads to a particularly good fit to the L and M flux. 
The visible image (lUeta et al.ll2000l ) shows an elongated elliptical nebula around the central 
star and suggests that one is looking at a bipolar nebula at some intermediate orientation. 

IRAS 06530—0213. — The SED of this object can be fitted by a normal, detached and 
extended cool dust shell with no extra hot dust component. The emission platea u from 8 to 



14 iim . and the strong 15.8 /xm emission feature were not fitted. The visible image (jUeta et al. 



2OOOI ) suggests that one is looking at a bipolar nebula at some intermediate orientation. 



IRAS 07430+1115. — The modeling of the SED of this object was more difficult than 
most. An adjustment of the density exponent to its best value resulted in a = 1, implying 
a decreasing mass loss rate. This is both unlikely and contrary to the results for the others, 
so we chose to fix a = 2, implying a constant mass loss rate. The resulting model fit was 
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only slightly inferior to that with a = 1. The SED was fitted with a low density envelope 
and a high density cool shell region. The observational excess in the L band suggests the 
presence of a hot component for this object also, but without an M-band observation, it 
is too unconstrained to model. The model fit is not as good as most of the other sources. 
The visible and near-infrared photometry cannot be fit by a normal model atmosphere with 
reddening; the resulting compromise fit includes more uncertainty in Ay than in most of 
the models. The model fit falls below the observations at 60 /im, is not good on the long 
wavelength side of the 21 /im features, and falls below the L-band fiux, as mentioned above. 
Again, some excess emission is observed from 30 to 36 fim and the emission pla teau from 



8 to 14 fim is not included in the model fit. The visible image (lUeta et al.ll2000l ) shows an 
approximately circular nebula around the central star. 

IRAS 22574+6609. — This source is very faint in the visible (V=21.3) and it lacks a 
visible spectrum, so we assumed a typical PPN stellar temperature of 5500 K. The model 
requires a much higher optical depth (tq) for this object than the others, by a factor of 
10—100. The fit of the model to the SED of this object is good from 18 to 28 /xm and 
then from 60 to 100 /xm. The emission plateau rises above the continuum from 8 to 18 /xm. 
Beyond 28 /zm, the IRS spectrum continues to rise. An attempt to fit the IRS spectrum 
to 34 /im leads to a very poor fit to the 60 and 100 /im data. Without spectra extending 
to 45 /im, it is hard to know exactly where the continuum lies. To fit the K, L, and M 
band data, a hot component was added to the model, which contributes less than 2% of 
the optical depth at 11.2 /im. While this brings these data into agreement, we are still left 
with poor agreement at the shorter wavelengths, with extra fiux observed below 2.2 /im 
and through the visible^! Such behavior in dust shell models has been seen previously in 
one-dimensional models of extreme carbon stars including IRC+10216 and AFGL 3068, and 
it is attributed to a non-spherical dust shell geometry. This may also be the cause in IRAS 



(Ueta et al. 


2000; 


Su et al. 


2001) 



due a low optical depth path (hole) out of the dus t shell through which some scattered 



escaping, similar to the situation in IRC+10216 (IMen'shchikov. Hofmann. &: Weigelt 



ight i s 



2002). 



This m ay be consistent with the off-center intensity peak seen by ISanchez Contreras et al. 



(120061 ) in their high resolution near-IR images of IRAS 22574+6609. They show that this 
peak is not at the center of the nebula and thus is not the central star seen directly. We have 
not tried to add extra complexity to the model simply to fit the visible wavelength portion 



^The near-infrared measurements of IRAS 22574+6609 by iHrivnak fc KwoM (|l991br ) made with a 15" 
aperture likely include the nearby (5" south) source 2MASS 22591874+6625430, which is brighter in J but 
fainter in H and K. 
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of the SED. These new model parameters are significantly different from our previous model 
for this source. However, those were based on much poorer quality ISO spectra, in which we 
chose to disregard the observed ISO spectrum from 27 to 45 /im ( iHrivnak. Volk. &: Kwok 
2000h. The new mode l was best fit with Ay = 0.0; however, galactic extinction studies by 
Neckel fc Klard (Il980l ) suggest a value >2.5, so we fixed Ay = 3.0. The fit was slightly worse. 



but at a level that is negligible on Figure [61 

IRAS 23304+6I47. — The model gives a good fit to the SED apart from the emission 
plateau shortward of 18 fj,m. A hot dust component is required to fit the M-band excess. 
This hot dust con tributes less than 2% of the optical depth at 11.2 /xm. The visible image 
( ISahai et al.ll2007l ) is similar to that of IRAS 06530—0213, and also suggests a bipolar nebula 
at some intermediate orientation. This new model for IRAS 23304+6147 diffe rs slightly from 
the one that we previously published based in ISO spectra (I Volk et al.l 120021 ). due primarily 
to the adoption of a stellar model rather than a black body for the central star. 

Thus it can be seen that the model fits are generally good (excluding the unmodelled 
AIB emission plateau from 8 to 14 //m), except for IRAS 22574+6609 and perhaps IRAS 
07430+1115. Observations in the 5—10 /im range are needed to better constrain the con- 
tinuum on the short wavelength side of the AIB emission plateau. Although most of these 
PPNs are known to have a bipolar morphology, the assumption of spherical symmetry will 
have little effect on the fitting of the mid-infrared emission and the calculated mass-loss rate. 
This is the case because the dust shells are optically thin in the infrared where the dust is 
detected. The only exception to this might be IRAS 22574+6609, where the derived optical 
depth is much higher than for the other objects. In this case, the morphology might affect 
the mid-infrared emission and the resulting observed flux would then depend upon observing 
orientation. 

The cool dust temperatures are similar in these models, with values in the range 120—150 
K at the inner radius of the cool dust shell. Including an additional six 21 fim sources 
modelled earlier based on ISO spectra, t he full range of Th (c) is 140 to 210 K (disregarding 
the earlier model of IRAS 22574+6609: Iflrivnak. Volk. fc Kwo^boooh . The absence of C- 
rich PPNs with higher values of Td(c) arises partly from the absence of cooler, spectral type 
K cent ral stars with small er detached dust shells, even though there is no bias in detecting 
these ( IVolk fc Kwoklll989l ). The absence of lower values of Td(c) at the inner radius of the 
21 /im sources suggests that either the carrier is not excited at these lower dust temperatures 
or that it is destroyed due to the presence of the harder UV radiation field found when the 
central star has evolved to high temperatures. The value of Td(c) observed in a particular 
source depends sensitively on the rate of evolution of the central star as well as the expansion 
velocity of the detached shell. The temperatures for the hot dust are in the range ~500— 800 
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K. For IRAS 05113+1347, 05341+0852, and 23304+6147, the density of the cool dust is 
10,000-15,000 times that of the hot dust, but for IRAS 22574+6609, the ratio is only ~200, 
due to the relatively larger amount of hot dust. 

Many of the values in Table [5] are scaled to the distance (D), and the mass loss rate (M) 
is also scaled to the expans ion velocity (V ). Assuming a luminosity of 8300 L©, appropriate 



for a core mass of 0.63 Mq ( lBlockeiill995al ) and published expansion velocities, the distances 
and mass loss rates for each object were calculated. These are listed in Table El The mass 
loss rates range from a 10~^ Mq yr~^ to as high as a 6 x 10~^ Mq yr~^. These later values are 
very high and cannot be sustained very long in medium and low mass stars. The total mass 
of the cool dust shells (AM) range from 0.22 to 0.95 Mq. The total mass of the enhanced 
shell is a more robust value than the mass loss rate since the width of the shell is not tightly 
constrained. The duration of the intense mass loss was calculated based on the shell radii 
and the expansion velocity, and ranges from 160 to 820 years; again, these values, based on 
the width of the shell, are not tightly constrained. Thus the results of the model fits to the 
SEDs imply that these post-AGB objects lose most of their ejected mass in a short duration 
of time, 10^—10^ years. 

A kinematical age since the end of the end of the high mass-loss phase was calculated 
based on the inner radius of the cool dust shell and the expansion velocity; this ranges from 
1000 to 3000 yr for five of the objects but is only 300 yr for IR AS 2257 4 +6609 . These values 



are consistent with the evolutionary ages for post-AGB stars iBlockerl (Il995bl ). The results 
suggest that IRAS 22574+6609 may be more massive than the other sources, based on its 
larger value of the ejected mass, the larger expansion velocity, and the apparently shorter 
time scale of evolution. The kinematical ages also provide a time scale for the formation of 
the carrier of the 21 fim feature. 

These new results for the mass loss rate and mass of the ejected shell can be compared 
with other published values determined in other ways. Here we have confined ourselves to 
two different approaches that also involve either discrete denser inner shells or steep density 
laws. 

Ueta and Meixner have recently modeled the mid-infrared emission from the dusty 
envelopes of two carbon-rich PPNs based on high-resolution images. They used a model 
that had an extended AGB envelope combined with an equatorially-enhanced discrete shell 
ejected during a later superwind phase. Both of these are sources with 21 fim emission, 
although they are among the brighter ones and were not included in our present study. For 
IRAS 22272+543 5, the values determined were M = 4 x lO"*' Mq yr^^ and At = 1500 years 



flUeta et al.l 120011 ). For IRAS 07134+1005 (HD 56126), for which their model was based on 



CO interferometric observations in addition to mid-infrared images, they determined M 



3 X 10 ^ Mq yr ^ and At = 840 years (iMeixner et al.ll2004j ). These mass loss rate values 
are a factor of 10 to 100 lower than what we have determined here. 



Hrivnak & Beiging (2005) determined mass loss rates and shell masses based on their 
measurements of CO gas in several transitional states, which allowed them to sample the 
CO emission further out and closer in to the star. The observations were modelled with a 
one-dimensional statistical equilibrium, radiative transfer code. They did not use a discrete 
shell with inner and outer boundary, as we have done here, but they did use a steep density 
law (a=3), which produces a dense inner shell. Their values for the mass loss range from a 
few X 10~^ Mq yr~^ to a few x 10~^ yr~^, a factor of 10 less than we find here. For 
three objects in common, our values for the mass loss rate larger by factors of two to 100. 

A comparison of our new model results for mass loss rates with two different studies 
reveals that our rates are a few to 100 times larger, but are confined to a thinner shell. As 
mentioned above, our mass loss rates are not as certain as we would like because the shell 
widths are not well constrained and assumed to be thin. Similar model fits could likely be 
produced with somewhat thicker shells with approximately the same total mass. Future 
higher resolution (0.1") imaging than is presently available in the mid-infrared will better 
constrain these shell widths. It is the total mass that is the better determined parameter 
from our models. Indeed, comparing the dust masses determined in these studies with our 
values of 0.7—2.8 x 10"^ Mq (AM/330 from Table [6]), one finds much better agreement. 
Meixner et al. determined a similar value of 0.8 x 10~^ Mq for IRAS 07134-1-1005, while 
Ueta et al. determined a value 20 times smaller, 0.03 x lO'^ Mq for IRAS 22272-^5435. If we 
compare with the Hrivnak & Bieging mass loss values based on CO gas and adjust with our 
gas-to-dust ratio of 330, we find dust masses 2—4 times smaller than ours, 0.24 x 10^^ Mq 
and 0.42 x 10"^ Mq for IRAS 07134+1005 and IRAS 22272+5435, respectively. However, 
if our gas-to-dust mass ratio is overestimated by a factor of 2—4, then the agreement would 
be excellent. Note that Meixner et al. determined a gas-to-dust mass ratio of 75 for IRAS 
07134+1005. 



6. RESULTS AND CONCLUSIONS 



Our targets were selected with the goal of completing the sample of carbon-rich PPNs 
studied with mid-infrared spectroscopy. However, th ree additional carbon-rich PPNs were 
identified in the past few years that were not included (iReyniers et al.ll2004l . 120071 ). Therefore 
the sample of carbon-rich PPNs that have been observed in mid-infrared spectroscopy is 
complete except for these three objects. Below are listed results based primarily on the 
targets of this study, supplemented by more general results for the sample. 
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1. Three new 21 /iin sources were discovered: IRAS 06530—0213, 07 430+1115 and 



1947 7+2401. This brings to 16 the total of 21 yum sources among the PPNs (IHrivnak et al. 



20081 ). Improved spectra of this feature were obtained for four others. 



2. The shape and central wavelength (20.1+0.1 yum) of the the 21 /xm feature agrees 
with that found previously for the brighter sources using ISO spectra. 

3. These spectra do not show the resolution of the 30 /xm feature into two features as 
was claimed based on ISO spectra. The suggested separation apparently was caused by a 
bad detector band in the ISO spectra and is not real. 

4. The 11.3 /im AIB feature is seen in all of the sources. The feature has a relatively 
consistent peak wavelength of 11.33+0.03 //m. However, the presence of the features at 12.3 
and 13.3 /im and their peak wavelengths and shapes vary among the sources. 

5. These results strengthen the correlation between the presence of the 21 /um feature 
and evidence of carbon- rich nature of the sources. In particular, all of the 21 /im sources are 
carbon-rich (' C/O>1.0), possess rn olecular carbon absorption features (C2, C3), and possess 



AIB features (IHrivnak et all 120081 1. 



6. All carbon-rich PPNs that have been observed possess the 21 /xm, 30 /xm, and AIB 
emission features. An exception may be IRAS 01005+7910 (OBe), the PPN with the hottest 
carbon-rich central star, which may or may not possess a weak 21 /xm feature. We also note 
that IRAS 19500—1709 is unusual in possessing the 2 1 /im (weak) and 30 /xm (strong) 



features but weak or few AIB features (jVolk et al.ll2002l ). 



7. The 21 /xm feature is present in all (but perhaps one) of the Galactic carbon-rich 
PPNs but is absent or weak in extreme carbon stars and PNs. This suggests that the 
carrier is a common part of the outflow during the PPN phase and that the radiation field 
is important to the presence of the feature. We propose (or This supports the idea) that 
the 21 /im carrier is produced in the extreme mass loss near the end of the AGB phase, is 
excited during the PPN phase, and is then destroyed in the more extreme UV radiation field 
of the subsequent PN phase. 

8. A narrow feature at 13.7 /im due to C2H2 is observed in the spectra of five objects; it 
is seen in one source in absorption, in three in emission, and possibly one in both emission 
and absorption. These represent the first observations of C2H2 emission in a post-AGB 
object. C2H2 absorption is commonly seen as a strong feature in the preceding AGB stage of 
stellar evolution. While most may be incorporated into larger molecules in the evolution to 
PN, at least some is still present in some region around the star in the PPN phase. Studies of 
C2H2 provide a potential diagnostic of the physical conditions in the circumstellar envelope 
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and should be followed up with high-resolution spectroscopy. 

9. At ~15.8 fim a broad (~1.3 /xm), unidentified emission feature is seen in the Spitzer 
spectra of four sources. This feature also appears to be present in ISO spectra of several 
additional 21 /im sources. It is particularly strong in the two Spitzer sources with the 
relatively strongest 21 fim feature. This suggests a possible correlation between the two 
features which, if confirmed, could be a valuable aid in the identification of the 21 /im 
feature. 

10. 1-D modeling of the SEDs indicates a dense cool (120—150 K) shell of dust around 
the stars and also low density hot (500—800 K) dust closer to the star. The mass-loss rates 
when the shells detached are high, ranging from a 10"^ Mq yr'^ to a 6 x 10"^ Mq jt'^ but 
the highest rates were not sustained for long. These result in shell masses of 0.2 to 1.0 Mq, 
ejected over intervals of 160 to 820 years. 

11. The results of the modeling of the dust shells imply ages since the termination 
of extensive mass loss of 1000—3000 years for all of the sources except IRAS 22574+6609. 
These ages for the expansion of the shell are consistent with the evolutionary ages of low- 



mass, post-AGB central stars (lBlockerlll995bl ). as they should be. IRAS 22574+6609 appears 



to have only relatively recently left the AGB, ~300 years ago. 
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Table 1: Spitzer Observing Log 



IRAS ID 


LH 




SH 




Date 


Time^ 


Date 


Time*^ 


05113+1347 


2006 Mar 17 


3 X 14 


2004 Oct 03^ 


3 X 30 


05341+0852 


2007 Apr 16 


3 X 14 


2004 Mar 22'' 


3 X 30 


06530-0213 


2006 Apr 25 


3 X 14 


2006 Apr 25 


3 X 30 


07430+1115 


2007 Apr 27 


3 X 14 


2007 Apr 27 


3 X 30 


19477+2401 


2005 Oct 16 


3x6 


2004 May 13*^ 


3x6 


22574+6609 


2005 Oct 12 


3x6 


2004 Aug 10^ 


4 X 30 


23304+6147 


2005 Dec 15 


3x6 


2005 Dec 15 


3 X 30 



"Number of scan cycles times the on-source integration time in seconds. 
''Observed in GT093, PI: D. Cruiksiiank. 
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Table 2: Target Positions (2000.0) 



IRAS ID 



2MASS ID 



RA 



Decl. 



05113+1347 
05341+0852 
06530-0230 
07430+1115 
19477+2401 
22574+6609 
23304+6147 



05140775+1350282 
05365506+0854087 
06553181-0217283 
07455139+1108196 
19495491+2408532 
22591835+6625482 
23324479+6203491 



05:14:07.76 
05:36:55.06 
06:55:31.82 
07:45:51.39 
19:49:54.91 
22:59:18.36 
23:32:44.79 



+ 13:50:28.3 
+08:54:08.7 
-02:17:28.3 
+11:08:19.6 
+24:08:53.3 
+66:25:48.2 
+62:03:49.1 
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Table 3: 21 //m Feature Profile and Relative Strengths 





Relative 








Feature 


Width (FWZL)^ 


\ b 


IRAS ID 


Strength 


(/xm) 


{/Jim) 


05113+1347 


0.49 


5.4 


20.2 


05341+0852 


0.35 


5.4 


20.1 


06530-0230 


1.12 


6.0 


20.1 


07430+1115 


0.21 


4.6 


20.2 


22574+6609 


0.37 


5.7 


20.1 


23304+6147 


1.46 


5.7 


20.1 



"Full width at zero level. 

''Wavelength of feature peak in a plot of F^, vs. A . 



Table 4. Summary of Spectral Emission Features 



IRAS ID 


SpT 




AIB Features (3-9 ;im)'^''= 


RcP^ 


11.3 fim 


21 


30 fira 


13.7 ^lm'' 


Other Features in Spitzer Spectra*^ 


05113+1347 


G8 la 


2.4 


3.3:, 3.4: 


1 


Y 


Y 


Y 


N 


12.1, 13.3: 15.8(mw):, 22.3 


05341+0852 


G2 0-Ia 


1.6 


3.3, 3.4, 6.2, 6.9, 8(br) 


1,2 


Y 


Y 


Y 


E 


12.4, 15.8(mw), (no LH Spitzer spectrum) 


06530-0213 


F5 I 


2.8 


3.3 


1 


Y 


Y 


Y 


E 


12.3, 13.2, 15.8(mw), 22.3 


07430+1115 


G5 0-Ia 




3.3, 3.4 


1 


Y 


Y 


Y 


E: 


12(mw) 


19477+2401 


... s 


... g 


f 


1 


h 


Y 


Y 


h 


(no SH spectrum), 22.3, 28.6:, 29.8:, 32.6: 


22574+6609 


... g 


... g 


■ • ■ f , 6.2, 6.9, 7.7, 8(br) 


1,2 


Y 


Y 


Y 


A 


12.6, 17(mw) 


23304+6147 


G2 la 


2.8 


3.3, 3.4:, 6.2, 6.9, 7.7, 8.6 (8br) 


1,3 


Y 


Y 


Y 


A/E: 


12.3, 13.2 15.9(mw), 22.3 



''References to C/O abundances: Multiple objects - I Van Winckel fc Revnier"3 l|2000l ') . Illeddv et alj l|2002h : IRAS 06530—0213 - iRevniers et al ] ll2004h . 
''Colons indicate an uncertain detection; mw = medium-wide feature; abs = absorption feature. 
'^Included also are aliphatic infrared bands at 3.4 and 6.9 /im. 

■^References for AIB features: ( 1) iHrivnak. Geballe. fc Kwokl | |2007|) ; (2) iHrivnak. Volk. fc Kwokl hOOd ): (3) review of ISO archival spectrum. 

"E = emission; A = absorption feature; N = no feature. 

^Not observed in the 3 /xm spectral region. 

^Faint, with no visible spectra. 

''Not observed in this spectral region with Spitzer. 



Table 5. Results of 1-D Model Fits to the Spectral Energy Distributions 



Cool Dust Hot Dust 



IRAS ID 


To 

(11.22 fim) 


To 

(V) 


a 


T, 
(K) 


(L0/kpc2) 


(K) 


rin{c)/D 
(10--^ pc/kpc) 


rout(c)/D 
((10-3 pc/kpc)) 


M/VD 

^ km s~-'-kpc~-'- ' 


Td(h)'' 
(K) 


nn{h)/D 

(10^3 pc/kpc) 


M/VD 

^ km s~ -'^ kpc~ ^ 


p{c)/p{h)^ 


(ISM) 


05113+1347 


0.0085 


1.2 


2.0 


5500 


170 


132 


4.2 


5.4 


4.5 X 10"'' 


714 


0.043 


3.4 X 10-1° 


12,500 


2.8 


05341+0852 


0.0067 


0.94 


3.5 


6500 


126 


146 


2.5 




1.0 X 10-8 


486 


0.117 


6.7 X 10-3 


15,000 


3.7 


06530-0213 


0.0052 


0.73 


3.0 


7000 


280 


130 


5.4 




1.6 X 10-« 










5.0 


07430+1115 


0.018 


2.5 


2.0^^ 


5500 


187 


139 


2.8 


3.3 


1.1 X 10-6 










1.5 


22574+6609 


0.39 


55 


2.0 


5500 


145 


148 


1.0 


1.6 


3.0 X 10-5 


765 


0.029 


4.9 X 10-3 


7700 


3.0<= 


23304+6147 


0.030 


4.2 


2.5 


5500 


320 


122 


5.0 


6.4 


2.0 X 10-5 


609 


0.083 


1.7 X 10-3 


10,000 


1.8 



'^Dust temperature at rin. 
'^Density ratio at rin(c). 
Fixed value. 
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Table 6: Model Results for Mass Loss Rates, Shell Masses, Ejection Durations, and Kine- 
matic al Ages 





pa 




M 




rin(c) 


AM 


t(AM)^ 


tkin 


IRAS ID 


(kpc) 


(km s~^) 


(Mo yr-i 


) 


(pc) 


(Mo) 


(yr) 


(yr) 


05113+1347 


7.0 


10 


3.2 X 10" 


-4 


0.030 


0.25 


820 


2900 


05341+0852 


8.2 


12 


1.0 X 10" 


■4 


0.020 


0.27 




1600 


06530-0230 


5.4 


14 


1.2 X 10- 


-4 


0.030 


0.22 




2100 


07430+1115 


6.7 


15 


1.1 X 10- 


-3 


0.019 


0.27 


180 


1200 


22574+6609 


7.6 


26 


5.8 X 10- 


-3 


0.008 


0.95 


160 


280 


23304+6147 


5.1 


12 


1.2 X 10- 


-3 


0.025 


0.66 


580 


2000 



"Distance based upon assumed luminosity of 8300 L©. 

''Expansion velocities as found in the literature and averaged in some cases. 
'^Duration of enhanced mass loss, determined from (rout(c) - rin(c))/V. 
'^Kinematical age, determined from rin(c)/V. 
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05113+1347 






10 15 20 

X\f X*J 


25 


30 35 


'7 05341+0852 


1 1 1 





10 15 20 


25 


30 35 


06530-0213 




, 1 , , , , 1 , , , - 


10 15 20 


25 


30 35 


07430+1115 






10 15 20 25 30 35 
Wavelength (|im) 



90 
60 
30 


60 
40 
20 


90 
60 
30 






10 15 20 


25 


30 


35 


22574+6609 


1 1 1 


1 1 1 












1 


10 15 20 


25 


30 


35 




15 20 25 30 35 
Wavelength (|im) 



Fig. 1. — Spitzer IRS high- resolution spectra of the seven sources observed with LH. They 
each display the 21 /im emission feature. 
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20 




9 12 15 18 21 

Wavelength (|im) 

Fig. 2. — Spitzer IRS high- resolution SH spectra of the six sources observed in this mode. 
They all show AIB emission features. The feature at 13.7 /im (attributed to C2H2) seen in 
three of the spectra is indicated; it may also be present in IRAS 23304+6147 and 07430+1115. 
(The noisy 14.5-14.6 /um spectral region has been removed from the spectra). 
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Fig. 3.— Comparison of Spitzer, ISO, and KAO spectra of IRAS 22574+6609 and IRAS 
23304+6147. 
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Fig. 4. — Comparison of the normalized 21 fim spectrum. (Excluded is IRAS 19477+2401 
due to the absence of the spectrum shortward of 18 fim with which to set the continuum.) 
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Fig. 5. — Comparison of the normalized 10—17 fim spectrum for the sources in this study. 
Included for comparison is the ISO spectrum of the young PN IRAS 21282+5050. Seen in 
all the spectra is the strong emission feature at 11.3 /zm; emission at 12.3 /im is seen in most 
and 13.3 and 15.8 /im in some of them. (The spectra are noisy in the range of 14.5 to 14.6 
fj,m.) 
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Fig. 6. — Fits of the models to the SEDs. Included are the ground-based (2MASS, DENIS, 
NIR=near-IR, MIR=mid-IR, Other), MSX, IRAS, and ISO data, in addition to the Spitzer 
IRS spectra. 



